Anderson transition in low-dimensional disordered systems driven by nonrandom 

long-range hopping 
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The single-parameter scaling hypothesis predicts the absence of delocalized states for noninter- 
acting quasiparticles in low-dimensional disordered systems. We show analytically and numerically 
that extended states may occur in the one- and two-dimensional Anderson model with a nonran- 
dom hopping falling off as some power of the distance between sites. The different size scaling of 
the bare level spacing and the renormalized magnitude of the disorder seen by the quasiparticles 
finally results in the derealization of states at one of the band edges of the quasiparticle energy 
spectrum. The delocalized nature of those eigenstates is investigated by numerical diagonalization 
of the Hamiltonian and by the supersymmetric method for disorder averaging, combined with a 
renormalization group analysis. 
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Localization of noninteracting quasiparticles in ran- 
dom media with time-reversal symmetry and finite-range 
hopping have been extensively studied since the semi- 
nal paper by Anderson Q. The hypothesis of single- 
parameter scaling, introduced in Ref. ||, led to the general 
belief that all eigenstates of noninteracting quasiparticles 
were exponentially localized in one (ID) and two (2D) 
dimensions (see Refs. j|| [| for a comprehensive review) 
and that localization-delocalization transitions no longer 
exist in the thermodynamic limit. Despite models with 
finite-range hopping work nicely in describing a variety 
of materials, long-range hopping are often found in dif- 
ferent physical systems (e.g. Frenkel excitons). Random 
long-range hopping was found to give rise to derealiza- 
tion of states not only in three dimensional systems 0] 
but in any dimension (|, ^, 0, Recent studies Q re- 
vised the validity of the single-parameter scaling hypoth- 
esis even within the original ID Anderson model with 
nearest-neighbor coupling, although did not question the 
statement that all eigenstates in ID random systems are 
localized. 

In this Letter we present analytical and numerical 
proofs that a localization-delocalization transition may 
occur in ID and 2D systems with diagonal disorder and 
nonrandom intersite coupling which falls off according 
to a power-like law. Apart from the importance of this 
finding from a general point of view, it may be relevant 
for several physical systems. As an example, let us men- 
tion dipolar Frenkel excitons on 2D regular lattices where 
molecules are subjected to randomness due to a disor- 
dered environment pQ |. Biological light-harvesting an- 
tenna systems represent a realization of the model we are 
dealing with [|ll} [l2| . Magnons in ID and 2D disordered 
spin systems provide one more example of interest. 



We consider the Anderson Hamiltonian on a d-dimen- 
sional (d — 1, 2) simple lattice with Af — N d sites: 



W = ^£ n |n)(n|+^ J nm |n)(m| 



(1) 



where |n) is the ket- vector of the state localized at site 
n, and {e n } are random site energies, assumed to be un- 
corrected for different sites and distributed uniformly 
within an interval [—A/2, A/2], thus having zero mean 
and standard deviation a — A/^/12. The hopping inte- 
grals between lattice sites m and n will be taken in the 
form J mn = J/|m — n| M (J mm = 0), where J > with- 
out loss of generality and the lattice constant is set to 
unity. We stress that hopping integrals do not fluctuate. 

For our qualitative reasoning we rewrite the Hamil- 
tonian ([j]) in the Bloch wave representation, |k) = 
A/" -1 / 2 J2 n exp(ikn)|n), with periodic boundary condi- 
tions. It then reads 



w = 5> k |k><k| + £(5wv|k><k'| 



kk' 



ri-ZO 



(<5W)kk< = ^£en^ k - k > 



(2a) 



(2b) 



(2c) 



where the wavenumbers k and k' run over the first Bril- 
louin zone. Notice that \i > d to ensure the convergence 



of (2b) in the thermodynamic limit. 

The key point of our qualitative arguments is as fol- 
lows EJ]. We compare the size scaling of the typical 



2 



magnitude of the scattering matrix (STi)uk' with the size 
scaling of the level spacing SE in the bare quasiparticle 
spectrum of the homogeneous Hamiltonian (A = 0). In 
particular, we focus our attention on those eigenstates 
laying close to the band edges k = (top) and k = 7r 
(bottom), where ir = tt and (n, 7r) for ID and 2D systems, 
respectively. The typical fluctuation of the scattering ma- 
trix ((JTi)kk' is °cff = o/N d / 2 . Thus, in spite of the fact 
that the magnitude of the disorder is a, the quasiparti- 
cle sees an effectively reduced value a c ft- It is important 
that CT e ff scales inversely proportional to N d / 2 . Straight- 
forward calculations of the bare energy spectrum ( pb| ) 
close to the band edges give the following results: 



£ k ~ E - JA d {n) |k|^- d , |k| -» 



JB d {y) |k - 7V\ 



(3a) 



(3b) 



where Eq and E^ are the band-edge energies and Ad{n) 
and Bd(fj,) are dimensionless constants. From @ it fol- 
lows that the level spacing scales as 6E ~ N~^ +d at the 
top of the band, while at the bottom one gets SE ~ N~ 2 . 

The matrix (<57i)kk' couples the bare (extended) quasi- 
particle states to each other and may result in their lo- 
calization within a region of size smaller than the system 
size. It seems reasonable to assume that the states will 
be weakly coupled and consequently will be delocalized 
over the whole system provided the inequality cr e g <C SE 
holds. It is remarkable that for fj, < id/2, the level spac- 
ing SE at the top of the band diminishes upon increasing 
N slower than the effective magnitude of disorder o~ e g. 
Therefore, if the coupling between bare states is weak 
for some finite N (cr c ff <C SE) then it will become even 
weaker upon increasing N. Consequently, one may ex- 
pect that the state will remain extended in the thermody- 
namic limit N — > oo . It is also reasonable to assume that 
disorder of magnitude larger than the bare bandwidth 
will localize all the states. From the above arguments 
we conjectured the existence of an Anderson transition 
in ID and 2D systems with diagonal disorder and non- 
random long-range hopping as long as /i < 3rf/2. Below 
we provide analytical and numerical confirmations of this 
conjecture. 

Concerning the parabolic range of the energy spectrum 
(close to the bottom of the band), we notice that the level 
spacing diminishes as A^ -2 upon increasing the lattice 
size, i.e., faster than the effective magnitude of disorder 
cr c ff. Now, even if o~ c ff <C SE for a small lattice size and 
the states are delocalized, the above inequality will be 
reverted for larger N, resulting finally in the localization 
of those states. The same conclusion holds for both band 
edges within the nearest-neighbor approximation, where 
the level spacing is always SE ~ N~ 2 . 

A supersymmetric method for disorder averaging |T4| , 
|l5[ , combined with a renormalization group (RG) anal- 
ysis, provide support to the above arguments. In short 



(the details will be published elsewhere), the sequence of 
our steps is as follows. As a first step, we consider the 
one particle Green's function with the fermionic partition 
function Zq and the bare action So defined as 



Z Q = J Y[ # n e~ So = det (H - SI) 



S = 1 V'n (n - £I)r 



Ipr, 



(4a) 



(4b) 



where £ = E + i0 + and I is the identity matrix. In- 
troducing bosonic ghosts (3,(3 and expressing 1/Zq as a 
path integral p5| , we then average the one particle Green 
function using the Gaussian probability distribution of 
site energies, P(s n ) — (1/ng) 1 / 2 exp(— e 2 /.g), instead of 
the box distribution introduced in the beginning. This 
allows us to perform the integration over site energies 
explicitly. The effective action 

Scff = i ^ [y>n(Jnm ~ £S nm )lp m 
n m 



(3n(Jnm £^nm)/^m] 



(5) 



which appears after averaging, will be the main object 
of our RG analysis. For doing this, it is convenient to 
rewrite the action in the k-representation and to regroup 
the terms as follows S e s = Skin + Ss + S g , where 



Skin = - iJA d J d d k Ikl^^k) ^(k) 

+ /5(k) /3(k)] , (6a) 



Sf 



i£ / d d k [ip(k) ^(k) + /3(k) (3(k)} , (6b) 



S n = 



r 4 

f- / TTd d k i <5(k 1 +k 2 -k3-k4) 

x^kx) ^(k 2 ) + ^(k 1 )/3(k 2 )] 

X[^(k3)^(k4)+^(k3)/3(k 4 )] . (6C) 



We have absorbed the constant Eq into E in d6b|) . Note 
as well that the integration over momenta are restricted 
to a d dimensional sphere of radius A, i.e. |k| < A, with 
A being an ultraviolet cutoff. 

The action @ is the starting point of our RG analysis, 
which is inspired by Shankar's approach |jl6f to fermionic 
condensed matter systems. The key observation is that 
the kinetic part of the action, Skin, is invariant under 
the following scaling transformation of the cutoff A, the 
momenta k and the fields: 



A ^ A' = A/6, 



k' = b k, b > 1 (7a) 
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(7b) 



For generic values of ji > d, Eq. ([7a]) is a non stan- 
dard scaling law which emerges from the unusual kinetic 
term (|6aj). 

Driven by (0), the mass term Ss also transforms into 
itself with a new coupling constant £' given by £' — 
b^~ d £. Thus £ is a relevant perturbation of the free 
action Skin, as it is always the case of mass terms |l6fl . 

The term S g also transforms onto itself under the fol- 
lowing RG transformation g' — b 2 ^ 3d g. This equation 
implies that g eventually goes to zero upon increasing b 
provided that /i < 3d/2. Hence, randomness vanishes 
in the low energy effective theory for fj, < 3d/2. On the 
contrary, g runs to stronger coupling whenever [i > 3d/ 2. 

Finally, the coupling g is marginal at tree level for fi = 
3d/2, and one has to consider the one loop effects to see 
its fate. This can be done using the techniques developed 
in Ref. Let us present our main results. The RG 

flows of the constants g and £, up to one loop, are given 
by 



ds 
d£ 
ds 



(2^-3d)g + f , 
(/i - d) £ - g , 



where s is the RG parameter defined as b 
g and £ are related to g and £ as follows 



9 = i-{JA d fK^ 

"d 



-3d 



(8a) 
(8b) 

: exp(s), and 
£ = ^JA d A^- d £, (9) 



with fi^ being the volume of the d-dimensional sphere. 

Equation ( |8a| ) has an unstable fixed point g* = 3d—2/i 
provided that /i < 3d/ 2 (see Fig. 0a). Below this point, 
the coupling g goes to zero asymptotically, while above g 
grows. For \i > 3d/2 the critical point disappears and the 
system always flows to strong coupling (see Fig. 0b). This 
signals about some changes in the density of states and 
the different nature of the eigenfunctions when passing 
from /j, < 3d/2 to /i > 3d/2. We would like to remark 
the fact that the critical value \x = 3c?/ '2 appeared in 
the present RG analysis coincides with that found on the 
basis of our qualitative arguments. 

„ g, = 3d - 2fJ 

a ) f < ¥ |—e= <£ <s • — => => 9 



b ) M > ¥ I 3> 3* 3* 3> 3* 5* 3* 3»- 3 

FIG. 1: RG-flows for the coupling constant g, Eqs. ( pa] ) 
and (|^), depending on the value of the exponent fi. 

It is usually argued that the averaged one-particle 
Green function does not carry information about the spa- 
tial extent of the eigenfunctions. Therefore, we chose 



the generalized inverse participation ratios (GIPR) as 
the relevant quantities to characterize localization prop- 
erties of the states. The standard definition reads I„ — 
(Sn l^nl 29 ), where "J^n is the probability amplitude of 
the normalized eigenstate v at site n, and the brackets 
denote disorder averaging. As is well known, the GIPR 
scales as 1^ ~ J\f~ q+1 for delocalized states (q ^ 1), 
while shows no size scaling for localized ones. 

Our main steps in calculating the GIPR are as fol- 
lows. We estimated the uppermost eigenfunction, \l/on, 
using perturbation theory with respect to the random 
term X)kk'(^) kk 'l' t )(k'| in (|2a|), namely considering the 
on-site energies small compared to the hopping parame- 
ter J, |e n | <C J. Then, \I>on can be written in real space 
as 



*0n 



(Ene 20 n) l/2 . 
m 

1 ^ e ik ( m ~ n ) 



(10a) 
(10b) 

(10c) 



In order to compute the GIPR we again made use of the 
the supersymmetric method for disorder averaging (as 
described above) with the bare action given by Eq. (||), 
as well as the replica trick introduced in Ref. [l7]. The 
latter reads 



rUl) 



lim — j 

r^O N d 



E 



,2</>„ 



(11) 



In doing so, we found that ~ N^^~^ d provided 
when d < [i < 3d/2; in other words, the generalized 
dimension equals the space dimension so that the upper- 
most state is delocalized, in full agreement with our qual- 
itative picture. In particular, notice that the so called 
inverse participation ratio (IPR) scales as 1^ ~ 

Since the previous analytical study of the GIPR was 
perturbative, we have also carried out a numerical study 
of the model to support the validity of our conclusions. 
We took advantage of the Lanczos method |l8| as well as 
the density matrix renormalization group approach [ jl9| , 
allowing one to calculate some few eigenstates of the 
Hamiltonian (|]) for rather large system size. In Fig. || 
we plotted the IPR (q = 2) of the uppermost state as a 
function of the system size Af = N d for different degrees 
of disorder A. The behaviour of the other top states is 
similar to that which manifests the uppermost state. Ob- 
serving Fig. |2| we conclude that the uppermost state is 
delocalized even for a moderately high value of the de- 
gree of disorder (A = 8J in ID systems and A = 40J 
in 2D system), provided d < [i < id/2. For comparison, 
the ID (2D) bandwidth for = 5/4 (jjt = 9/4) in the ab- 
sence of disorder is of the order of 10. 5 J (28 J). However, 
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for large degree of disorder the IPR remains constant 
on increasing the system size (see Fig. ^ for fx = 5/4, 
d = 1 and A = 40 J), indicating that the uppermost 
eigenstate is localized. Therefore, the top eigenstates un- 
dergo the Anderson transition on increasing A whenever 
d < fi < Zd/2. It is to be noticed the absence of scaling 
of the IPR and the subsequent localization for d = 1, 
fj, = 3, and A = 8J. This result is in full correspondence 
with the analytical analysis stating that no transition is 
expected in this case. 
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FIG. 2: Scaling of the IPR of the uppermost eigenstate as a 
function of the number of sites Af = N d for different magni- 
tudes of disorder A in ID (d = 1) and 2D (d = 2) lattices. 

In summary, we have shown analytically and numer- 
ically that a power-like nonrandom intersite hopping, 
Jnm = J/|n— m| M , may act towards derealization of 
quasiparticle states in low dimensional systems. In par- 
ticular, the states of the top of the band may be delo- 
calized at rather high magnitudes of disorder (A ^> J) 
and undergo a localization-delocalization transition as 
the magnitude of disorder increases. Scaling arguments 
provide a clear physical picture of the underlying mech- 
anism responsible for the derealization of the states, in 
spite of the low-dimensional (d < 2) geometry of the 
system. The different size scaling of the bare level spac- 
ing, SE ~ 7V _Al+d , and the magnitude of disorder seen 
by a quasiparticle, <r e g ~ A/N d / 2 , is the feature of the 
model responsible for this unusual behavior. SE de- 
creases slower than a e g upon increasing the system size 
as long as d < /i < 3d/ 2, resulting in the derealization 
of the corresponding quasiparticle states in the thermo- 
dynamic limit. We stress that the main finding of our 
study, namely the existence of the Anderson transition in 
a physically relevant model, has been concluded on the 
basis of three different approaches and the conclusions 
obtained are self-consistent. Most important, the validity 
of the scaling analysis is not limited to the present model. 



Indeed, it is established on solid grounds that the stan- 
dard, three dimensional Anderson model manifests the 
localization-delocalization transition at the band center. 
Within this model, the bare level spacing at the band 
center diminishes proportionally to A" 1 , while the mag- 
nitude of effective disorder goes down faster, ~ A~ 3 / 2 , 
thus being unable to localize the states at the band center 
for a moderate disorder. A strong disorder (large com- 
pared to the band width) localizes the states, giving rise 
a localization-delocalization transition. We can then be 
confident that this kind of scaling arguments may provide 
physical insight in several localization problems. 
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